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'. ABSTRACT 

We measure the distribution of velocity dispersions of groups of galaxies iden- 
tified in the UZC catalog; we use the distribution to derive the group mass 
function. We introduce a new method which makes efficient use of the entire 
magnitude limited catalog. Our determination of n(> o~t) includes a significant 
contribution from low luminosity systems that would be missing in a volume 
limited sample. We start from a model for the probability density function of 
the total number of group members and reproduce the observed distribution. We 
take several effects like the local fluctuations in volume density, limited sampling 
and group selection into account. We estimate the relation between total number 
of members, total luminosity and true velocity dispersion. We can then reproduce 
not only the observed distribution of a v but also the distributions of the number 
of group members, the total groups luminosity, and virial mass. The best fit to 
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the data in the true velocity dispersion range 100 kms < a T < 750 kms 

-3.4+ 1 - 3 

is a power law model with a slope n{> a T ) oc a T ' ~ 16 and a normalization 
n(a T > 750 kms' 1 ) = (1.27 ± 0.21) x 10' 5 h 3 Mpc~ 3 where a T is the true ve- 
locity dispersion of a group. The predictions of CDM models are consistent with 
the mass function we derive from this distribution of velocity dispersions. 

Subject headings: galaxies: clusters: groups: general — - cosmology: observations 
— cosmology: theory — - large scale structure of the universe 

1. Introduction 

Groups of galaxies contain approximately half of the galaxies within a magnitude limited 
redshift survey (see e. g. Ramella, Geller & Pisani 2002; Ramella, Pisani & Geller 1997; 
Trasarti-Battistoni 1998; Giuricin et al. 2000; Tucker et al. 2000; Adami & Mazure 
2002; Merchan & Zandivarez 2002). In hierarchical structure formation theories, these 
abundant systems are the natural progenitors of the galaxy clusters. Thus measurement 
of the distribution of the physical parameters of groups provides part of the fundamental 
basis for understanding the large-scale structure of the universe (Zabludoff & Geller 1994; 
Nolthenius, Klypin & Primack 1994; Diaferio et al. 1999; Carlberg et al. 2000; Girardi & 
Giuricin 2000). 

The velocity dispersion of a system of galaxies is an indicator of the depth of the potential 
associated with the system. Hence the distribution of velocity dispersions is the basis for 
the determination of the group mass function. Although the relation linking the velocity 
dispersion with the gravitational mass depends on the relative distribution of visible and dark 
mass and on the dynamical state of the galaxy system (Diaferio et al. 1993, but see also 
Xu, Fang & Wu 2000), the velocity dispersion itself is a physically interesting parameter. 
The data readily provide an estimate of the group velocity dispersion which is robust against 
the inclusion of faint members (Ramella, Geller & Huchra 1995; Ramella, Focardi & Geller 
1996; Mahdavi et al. 1999). The small median number of group members introduces a large 
scatter in the observational estimate, a v , of the dispersion of a system relative to the true 
value, or • This finite sampling effect also introduces distortion in the distribution of a v . 

Several investigators have determined n{> a T ) for optical galaxy clusters (e. g. Bahcall 
et al. 2002; Reiprich & Bohringer 2002; Borgani et al. 1999; Girardi et al. 1998; Fadda 
et al. 1996; Mazure et al. 1996; Bahcall & Cen 1993; Biviano et al. 1993). In addition, 
there are also measurements of the cluster x-ray temperature function (Ikebe et al. 2002; 
Pierpaoli, Scott & White 2001; Blanchard et al. 2000; Henry 2000; Markevitch 1998; Henry 
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& Arnaud 1991; Edge et al. 1990). On the theoretical side, large numerical simulations 
provide estimates of the predicted group/cluster mass function (e. g. Jenkins et al. 2001; 
Governato et al. 1999; Diaferio et al. 1999; Lacey & Cole 1994; Thomas et al. 1998). for 
a variety of cosmological models. 

For galaxy groups with low luminosity and hence low velocity dispersion, some deter- 
minations of n(> <jt) rest on small samples of systems in volume limited catalogs (e. g. 
Zabludoff et al. 1993). Others contain residual biases resulting from the magnitude limit of 
the samples (e. g. Moore et al. 1993; Girardi & Giuricin 2000; Heinamaki et al. 2003). 
In principle, restriction to a volume limited catalog mitigates the problem of selection effects 
introduced by the apparent magnitude limit of the redshift survey (Zabludoff et al. 1993). 
However, low luminosity systems are virtually absent from volume limited catalogs. 

The determination of n(> a?) from magnitude limited samples requires correction for 
the selection function. The standard approach is to weight each group by the maximum 
volume where a group could still be recognized as a triplet. The weighting scheme assumes 
that throughout the mass range the catalog covers a fair sample of the universe. These 
determinations also assume that the mass-to-light ratio, M/L, is independent of velocity 
dispersion. 

A reliable determination of the mass function of groups is particularly interesting be- 
cause there are still significant discrepancies between theory and observations at the low-mass 
end of the distribution (Martinez et al. 2002; Heinamaki et al. 2003). Standard observa- 
tional approaches to the determination of n{> a T ) are most uncertain at the low-mass end. 
For example Zabludoff et al. (1993) lack low-mass systems because of the absolute mag- 
nitude cut of their volume limited sample. For larger samples, the assumption of constant 
number of groups per unit volume is a potential problem. Even in these larger surveys, low 
mass systems probe a small volume where the distribution of systems is not homogeneous. 

Here we introduce a method for constructing n(> <jt) from all of the groups iden- 
tified in the magnitude limited UZC redshift survey. We improve the statistics of the 
n(> cry) determination in the range 100 km s" 1 < o T < 750 km s -1 . We determine 
n(> a T ) in a homogeneous way throughout the entire range in Ot ■ We account for the 
effects of limited sampling on the estimated velocity dispersion of each group and for local 
fluctuations in the number density of groups. 

We demonstrate that in the UZC, the intrinsic relations among the luminosity, velocity 
dispersion and richness of groups are well-represented by power laws. By using these model 
relations we can satisfactorily account for the observed distributions of all the group param- 
eters. These scaling relations for groups may be useful for evaluating the evolution of group 
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properties in catalogs extracted from deeper surveys. 

Section 3 discusses the group catalog. Section 4 outlines the main problems in esti- 
mating n{> <jt) and the basic method we apply. Section 5 discusses a power-law model for 
the relation between group richness and velocity dispersion. Section 6 describes the entire 
analysis method in detail and shows how we use the relation in Section 5. Section 7 con- 
tains the distribution of velocity dispersion and the group mass function. In Section 8 we 
discuss the scaling relations and demonstrate the broad consistency of the method we use. 
We summarize in Section 9. 



2. The group catalog 

We base our analysis on the UZC catalog (Falco et al. 1999) in the region defined by 
-2.5° < <5 1950 < 50° and 8 h < a 1950 < 17 h in the North Galactic Cap and by 20 h < a 1950 < 4 h 
in the South Galactic Cap. In this region the completeness in redshift is 98% . We discard 
the region —13° < b < 13° because of the greater Galactic absorption there (Padmanabhan, 
Tegmark & Hamilton 2001). 

Ramella, Geller & Pisani (2002) describe the catalog of UZC groups and list the mem- 
bership of all of the groups analyzed here. Ramella, Geller & Pisani (2002) identify the UZC 
groups with a friends-of-friends algorithm (FOFA) (Ramella, Pisani & Geller 1997; Huchra 
& Geller 1982). This algorithm is currently in wide use (Adami & Mazure 2002; Merchan 
& Zandivarez 2002; Girardi & Giuricin 2000; Tucker et al. 2000; Trasarti-Battistoni 1998) 
along with some applications of the dendogram analysis applied by Materne (1979) and by 
Tully (1987). The groups Ramella, Geller & Pisani (2002) identify are number density 
enhancements with 5pn/pn > 80 in redshift space (a, 5, V). 

The FOFA associates galaxies with a projected separation on the sky less than D = 
0.25/i _1 Mpc and a line-of-sight velocity difference less than U = 350 kms" 1 at a reference 
fiducial velocity v F = 1000 km s -1 . We scale the linking parameters D L = R(V) x D 
and Vl = R(V) x V with a function R(V) according to the prescription of Huchra & 
Geller (1982). We assume that the UZC galaxy luminosity function (LF) is in the Schechter 
(1976) form and adopt M* = -19.1, a = -1.1, and 0* = 0.04 h 3 Mpc~ 3 . We obtain 
these parameters by convolving the LF determined by Marzke et al. (1994) for a very 
similar sample with a Gaussian 0.3 magnitude error. These parameters improve those used 
by Ramella, Pisani & Geller (1997) to identify groups within the CfA2N survey. The 
differences between the Ramella, Pisani & Geller (1997) groups and the groups we identify 
now within the same region are negligible. 
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Within the UZC we identify 301 groups with at least 5 members and with mean velocities 
SOOkms^ 1 < V < 12000kms _1 . At the minimum group radial velocity V m i n = 500 kms -1 , 
the apparent magnitude limit of UZC, mu m = 15.5, corresponds to an absolute magnitude 
Mu m = mn m — 25 — 5 \og(y min / H ) = —13.00 + 51og(/i) and, with h — 1, to a luminosity 
\og(L Hm /L G ) = 7.39. 

To eliminate 39 groups with very low intrinsic luminosity, we limit our analysis to 
groups with mean velocity V > 3000 km s -1 . We refer to the remaining Nq = 262 groups 
as the UZCGG sample. These groups have true velocity dispersions larger than ~ 2.5 times 
the errors in the individual redshift determinations. The lower limit in the true velocity 
dispersion of our 262 groups is comparable with the lower limit of the velocity dispersion of 
the x-ray emitting groups detected within UZCGG (Mahdavi et al. 1999). 

In a group catalog selected from a redshift survey with a FOFA, some fraction of the 
groups are accidental superpositions. We have two measures of the fraction of true physical 
systems in the UZCGG. Ramella, Pisani & Geller (1997) use geometric simulations of the 
large-scale structure in the northern UZC region to demonstrate that 80% of the groups 
are probably physical systems. Diaferio et al. (1999) compare ACDM simulations with 
the northern portion of the UZC and conclude that the linking parameters we choose for 
group selection are in the optimal range for reproducing the catalog derived from n-body 
simulations. In that range the fraction of real groups is 70-80%. Mahdavi et al. (2000) 
cross-correlate a large portion of the UZCGG with the ROSAT All-Sky Survey (RASS). 61 
groups in the Mahdavi et al. (1999) sample have associated extended x-ray emission. The 
presence of hot x-ray emitting gas bound in the group potential well is a confirmation of the 
physical reality of the system. Mahdavi et al. (2000) use the groups detected in the x-ray 
to show that a minimum fraction of 40% of the groups in the UZCGG subsample are similar 
x-ray emitting systems undetectable with the RASS; thus they set a lower limit of 40% on 
the fraction of real physical systems in the group catalog. At least 40% and probably 70-80% 
of the UZCGG systems are real and it is thus reasonable to use their properties to derive 
physical constraints on the nature of groups of galaxies. 



3. The velocity dispersion distribution: the method 

The group velocity dispersion is a fundamental quantity both for studying the internal 
dynamics of the group and for understanding the processes which produce these galaxy 
systems (Frenk et al. 1990; Moore et al. 1993; Weinberg & Cole 1993; Zabludoff et 
al. 1993; Crone & Geller 1995). It is therefore important to have a reliable estimate of 
the distribution of velocity dispersions for groups, the most common bound systems in the 
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universe. 

Our goal is measurement of both the probability density of groups with a given velocity 
dispersion, f(a T ), and the space density of groups as a function of velocity dispersion, 
n(> <Jt) ■ One direct approach to the problem is analysis of a volume limited catalog of 
groups. This procedure significantly reduces the number of groups in the sample and, more 
importantly, it discards a large number of low luminosity systems largely with low velocity 
dispersions. For example, within the UZCGG there are 34 groups within the volume limited 
sample with velocity limit V* = 8300 km s _1 , and only 20 groups within V M = 12000 km 
s" 1 . 

Undersampling the internal velocity field of each system also affects the determination 
of n(> a?) f° r groups. The UZCGG groups have a median of 7 members: the estimate of 
their velocity dispersion o v (corrected for radial velocity uncertainty according to Danese, De 
Zotti & di Tullio 1980) has a large scatter around the true underlying velocity dispersion, 
ot ■ We use the unbiased estimator of Ledermann (1984) to compute a v . To illustrate 
the scatter, we use a Monte Carlo simulation to extract 1000 groups at random from a 
specified true distribution of cr T . We simulate the "observed" group velocity dispersion a v 
under the assumption that peculiar velocities within the group are distributed according to 
a Gaussian with true dispersion ot- The distribution of mean radial velocity and richness of 
the simulated groups is the same as for the groups observed in the UZCGG sample. Figure 
1 shows the relation between the true ot and the Monte Carlo sampled a v (the straight line 
represents equality of the two measures). The spread is quite large; it decreases at larger a v 
because the number of members observed is typically larger. For poorly sampled systems, 
the observed value of a v is a poor estimator of the true ot- In conclusion, the estimate of 
the distribution of the velocity dispersion is severely affected by several observational and 
sampling effects. 

To avoid these difficulties we assume that the peculiar velocities of group member galax- 
ies are gaussian-distributed. This assumption is warranted by our deeper sampling of groups 
with associated x-ray emission (Mahdavi et al. 1999). With an average of 35 galaxies per 
group, the velocity distributions are consistent with a Gaussian. 

We simulate the observed value of a v by Monte Carlo-sampling of the Gaussian velocity 
field with dispersion ot • We use the data to demonstrate that there is a relation between 
a T and the absolute group richness N. We confirm the relation by examining a set of 43 
more deeply sampled groups from Mahdavi et al. (1999) and Mahdavi & Geller (2001). 
We account for the selection effects introduced by: a) the apparent magnitude limit of the 
galaxy catalog, b) the distance dependent criterion of inclusion of a group in our analysis: 
Nmem > N min = 5, and finally c) fluctuation in the space density of groups p(V) . 
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We select groups according to the number of members. We therefore must estimate the 
probability density function, p(N), of the group absolute richness N, i.e. the total number 
of group members brighter than Ln m (or M[ im = —13.00). We use a power law model for 
p(N) and constrain both its slope and normalization. From p(N) we predict the number 
density of groups with more than N members. 

By using the relation between group richness and velocity dispersion along with the 
estimate of p(N), we can compute both /(cry) and n(> ax) ■ From these two distributions 
we can use Monte Carlo simulations to predict the observed distribution of a v . We can also 
predict the distribution of observed group luminosities, L mem , provided that the luminosity 
function is universal. 

The procedure includes the following steps: 

step 1: we estimate the relation between the number of group members (richness) N and 
the true velocity dispersion ax (Section 5) 

step 2: we choose a model for the probability density function of the group richness N: 
p(N) (Section 6) 

step 3 we estimate the group selection function by using p(N) (step 2) and then, by compar- 
ison with the data, we obtain the radial distribution of group density p(V) (Section 6) 

step 4: by using the N vs. a T relation (step 1), the model for p(N) (step 2) and the estimate 
of p(V) (step 3), we compute the cumulative distribution of the velocity dispersions 
n(> a v ) and compare it with the real data in order to assess the agreement with p(N) 
from step 2 (Section 7) 

step 5: by using the relations among richness, total luminosity, ax and the group mass, 
we compare the model distributions of all these quantities with the observations (Sec- 
tion 8). 

The main advantages of this procedure are; a) we determine n(> ax) for groups with Lx 
an order of magnitude fainter than for a volume limited catalog of groups identified within 
the same region b) we need not assume a constant space density of groups, and c) we use a 
sample of groups an order of magnitude larger than the the volume limited sample, and d) 
the mass function estimate does not require the assumption of constant M/L. 
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4. Group richness and velocity dispersion 



Several investigators have explored the relationship between various measures of "rich- 
ness" and the velocity dispersion of rich clusters of galaxies (Bahcall & Cen 1993): richer 
clusters tend to have larger velocity dispersions. Bahcall (1988), for example, fits a power- 
law relation between between the galaxy count within 0.25 h _1 Mpc and the system velocity 
dispersion for 23 Abell clusters and a sample of groups (Turner & Gott 1976). More re- 
cently, Yee & Ellingson (2002) have shown that the richness of CNOC1 clusters of galaxies 
is well correlated with their main physical parameters as derived from both optical and x-ray 
observations. Here we use the UZCGG data along with a set of well-sampled x-ray emitting 
groups (Mahdavi et al. 1999; Mahdavi & Geller 2001) to examine the relationship between 
group richness and velocity dispersion. 

We expect some relation between the total number of group members N and the true 
velocity dispersion, <tt, because N and Lt are related through the luminosity function, Lt 
and mass are related for bound systems, and mass and velocity dispersion are also related (if 
light traces mass). We do not require the same M/L for groups at every velocity dispersion. 
As a simple approximation we choose a power-law model for the relation between n and gt'- 



where N and a are scale factors and a is the exponent. 

To estimate the parameters Gq and a, we fix an arbitrary normalization in number, 
N = 100, and proceed in the following steps: 

1. We select the i — th group from the sample at random. This group has N rnem (i) 
members at a mean radial velocity V(i) and a total number of galaxies brighter than 
Mu m = —13.0, N(i) = N mem (i)/v(V(i)). v(V(i)) is the fraction of galaxies detected 
in a group at mean radial velocity V{i): 



The luminosity L cut (V(i)) corresponds to the UZC apparent magnitude limit at the 
radial velocity V(i) and 4>(L) is the luminosity function. 

The true velocity dispersion of the group is, according to our model, 




(1) 





(2) 



f Ll 2<f>(x)dx 




(3) 
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2. We sample N mem (i) velocities from a Gaussian with standard deviation ar(i) and 
compute the simulated velocity dispersion a simjV (i). 

3. We repeat these steps until we obtain the desired number of simulated groups. 

4. We perform a KS test between the simulated distribution of <J S im,v and the distribution 
of observed a v 

5. We vary cr and a in the interval a <G [1; 3] and <7 ( kms -1 ) G [400; 600] in a 21 x 21 
grid, and maximize the KS significance level. 

We iterate the previous five steps ten times with 5000 groups in each iteration and 
compute the median of the KS-test significance Sks of the comparison between the predicted 
and the observed distributions. 

With this procedure applied to the UZCGG, we obtain a = 510 kms -1 and a = 
1.8 with median S KS (<J V ) = 0.44. The 99% confidence interval is (1.4; 2.1) for a and 
(460 kms -1 ; 580 kms -1 ) for a . The value of a is sensitive to the definition and scale of 
"richness." 

Figure 2a shows the contours of constant Sks(&v) in the parameter space. The thick 
outermost isopleth is the 99% confidence level, the middle contour traces the 90% confidence 
level, and the inner thin contour traces the 50% confidence level. Note that the values of (Jq 
are scaled for group richness estimated at V > 3000 kms" 1 and hence for members brighter 
than —16.89. The corresponding value of the best fit cr scaled to V > 500 kms" 1 and 
M Um = -13.00 is 230 kms" 1 . 

To test the model in equation eq. 1 further, we obtain o" and a for a set of X-ray emitting 
groups kindly provided to us by A. Mahdavi (1999). Mahdavi et al. (2000) identified a set 
of groups with extended x-ray emission by cross-correlating a portion of the UZC with the 
ROSAT All-Sky Survey. The extended x-ray emission essentially guarantees the reality of 
the system. Mahdavi (1999) and Mahdavi (2003) obtained deeper redshift surveys of the 
groups and increased the typical membership to 30 galaxies. We use these enhanced data 
here to test the scaling relations derived for the UZC groups. 

There are 43 groups in the Mahdavi sample with mean radial velocities V > 5000 km s" 1 
and galaxies brighter than m < 16.5. The average number of group members is 20 and hence 
the value of the velocity dispersion is estimated more reliably than for UZCGG. Moreover, 
the diffuse X-ray emission strongly reduces the false group identification problem. Figure 2b 
shows the KS-test confidence level contours for the X-ray sample. The best fit value for the 
slope a is 1.7. The 99% confidence intervals are (1.0, 2.7) for a and (410 kms -1 , 590 kms -1 ) 
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for a . Because of the different completeness limits and different cuts in mean radial velocity, 
the values of the o"o normalization in the two plots for X-ray and UZCGG differ by a factor 
1.06. This offset improves the overlap of the contours in the <7o direction. Although the area 
covered by the contours is quite large, it is quite remarkable that X-ray groups yield values 
for a and cr Q nearly coincident with the UZCGG sample. 

Figure 3 shows N as a function of a v for the x-ray sample. The solid line is the result 
we obtained with the Monte Carlo procedure. The thin line is the bisector of the two fits 
obtained from a straightforward x 2 fit to N vs a v and to a v vs N respectively. We find 
log(iV) = 1.27l 42 x l°g(°"f) ~~ 1-47^2 04 consistent with the Monte Carlo approach. The 
Monte Carlo yields a steeper slope because it weights less the rare groups with very low 
velocity dispersion which also depart from the L x — a scaling relations (Mahdavi & Geller 



Finally we note that our result agrees with the power law relation that Yee & Ellingson 
(2002) find between richness parameter, B cg , and radial velocity dispersion, a"i, for CNOC1 
clusters. Yee & Ellingson (2002) find B cg oc ai L8±a2 within the range 600 kms" 1 < a x < 
1300 km s -1 . The agreement is interesting also because there is only a marginal overlap 
between our high a range and the low a range of CNOC1 clusters. 



Here we investigate the probability density function of groups, p(N), for the UZCGG 
sample. Most (92%) of the groups have N < N r = 840 (or a T < a r = 750 kms^ 1 ). Sparse 
sampling of the richest systems (clusters) limits the range over which we can reliably estimate 
p(N) to 20 < iV < 840. In our determination of p(N), we include all groups. Groups with 
with N > N r are visible (i.e. N j, s > 5) throughout the whole sample volume. These groups 
have marginal weight in the determination of the slope of p(N) but they are important for 
the normalization of n(> (Jt) ■ 

We assume a power-law model for p(N) : 



where N norm is a scale factor and A is a normalization constant. 

Here we set N norm = N min , with N min = 5 (see the previous section). We note that the 
choice of N norm is arbitrary because the power-law is scale free. 

The Press & Schechter (1974) formalism, PS hereafter, provides a basis for this assump- 



2001). 



5. The probability density function of groups 




(4) 
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tion. Valageas & Schaeffer (1997) argue that in the low mass range occupied by poor galaxy- 
systems the probability density function p(N) is well described by a power law. Recently 
Jenkins et al. (2001) use numerical simulations to obtain similar results. We make the 
argument appropriate for our data here. 

The PS model describes the mass function as: 



n PS (M)dM = K 




where K is a normalization factor, M* is a scale factor and the exponents a and (3 are 
linked to the spectral index, n e //, of primordial mass fluctuations: a = n e ///6 — 3/2 and 
(3=1+ n e ff/3. If we suppose that the group mass M is a power law function of the group 
number of galaxies N, 

M oc N k (6) 
we can rewrite the mass function in terms of N: 



The local slope of the PS model is: 
where x = N/N*. In other words: 

7 = (3kx 0k -ak-k+1 (9) 



The slope 7 depends on x and hence on N, but for the range of interest here (20 < 
iV < 840), 7 goes from 2.2 to 3.3. This range is in good agreement with the 90% confidence 
interval in our fit, namely (2.7, 3.2) (see section 7). By taking iV as the average value < N > 
of the N estimated for the observed groups, we can conclude that Equation 9 links the two 
parameters n e ff and N+ that characterize the PS model. Given a value of n e ff we should 
obtain a "best fit" N+: 



8k \ 1/m 

N*=<N>[ (10) 

1 7 + ak + k- l) K ' 



is shown in Figure 4. 
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Over the range covered by our data, the Press-Schechter model is well-approximated by 
a simple power law model for the mass and N functions. 

We next evaluate 7 in Eq. 4 from the observed group catalog. We start by estimating 
the probability of finding a group at mean velocity V with at least five observed members, 
Nmern > N min = 5. We call this function the group selection function, £(V). 

The probability that a group has at least N total members is: 

r+oo 

S(N) = / p(x)dx (11) 

J N 

with S(N) normalized so that S(N norm = N min = 5) = 1. 
Consequently, the group selection function £(V) is: 



W) = S (^y) (12) 
Clearly, S(V) depends on both <p{L) and p(N). 

From the group selection function, we can compute the distributions of V and N mem 
from the model for p(N). The distribution of group radial velocities is: 



r v 

C(V) = / ACLv 2 p(v)X(v)dv (13) 
Jv m 

where V m = 3000 kms^ 1 , Ail = 3.16 sr is the solid angle of UZCGG, and p(V) is the 
density of groups as a function of their mean radial velocity. The distribution of observed 
members is: 

Q(N mem ) = A{lv 2 p(v)S (^) dv (14) 

with Vm = 12000 kms -1 , the redshift limit of the group catalog. 

We compare the simulated C(V) and Q(N mem ) with the observed distributions of the 
same quantities. To make the comparison we must solve a system of two equations in the 
two functions p(N) and p(V). We have an assumed form for p(N). To estimate p(V) we 
invert the equation for C(V) (Eq. 13). We approximate p(V) with a sequence of N s h e ii radial 
velocity shells, each with constant density, i.e.: 



p(V) = pi if K-i < v < v t 



(15) 



with % — 1, . . . , N s heii- If there are 5Ni observed groups within the % — th shell, pi is: 



(16) 



Pi 



Anj;;_ iV ^(v)dv 



We choose the limits of the shells, V$, so that all shells contain the same number of 
groups. We start with N s h e ii = 1 and perform a KS-test to estimate the significance Sks(V) 
of the agreement between the model prediction C(V) and the observed distribution of V. 
We increase N she u until S KS (V) > 0.5. 

We have thus derived a reasonable approximation to p(V) and we next use it in eq. 14 
to compute Q(N mem ). We use a KS-test to evaluate the agreement between Q(N mem ) and 
the observed distribution of N mem . If the agreement is satisfactory, we conclude that we 
have a satisfactory slope for p(N). If, on the contrary, the distributions are inconsistent, we 
repeat the whole procedure with a different slope for p(N). Once we determine p(V) and 
p(N) , we can compute the total abundance of seen and unseen groups richer than a fixed 
threshold within the UZCGG volume, i.e. the group multiplicity function p{> N). 

The total number Af G of groups within UZCGG is: 



and the group average density is: 



where V = 1.79 x 10 6 h 3 Mpc 3 is the UZCGG volume. Under the assumption that p(N) 
does not depend on position within the UZCGG volume, the multiplicity function is: 



For groups with N large enough to be observable throughout the UZCGG volume, we 
can compare our model prediction p(N) directly with the observations. We call these groups 
robust and use them in order to normalize the function p(> N). 

To go fromp(iV) and p(> N) to n(> a T ) we use the relation between N and <jt (equation 
1). The data constrain the power-law slope a and the normalization er . The relations we 
need are: 




(17) 
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f(a T )=p(N) 



dN 



(20) 



and 



r+oo 



n 



(> °t)=p\ f(x)d. 



X 



(21) 



We can compute n(> <jt) by inserting equation 4 into equation 20. 



Because ot is related to N as in the eq. 1, we can write 



\ "norm / 



(22) 



with 7 S = 7a — a + 1 and cr. 



norm 



is a normalization scale. 



To estimate the value of the exponent 7 we compare the observations with the predic- 
tions of the p(N) model obtained by inserting eq. 4 for p(N) into eq.s 13 and 14. The best 
fit value is 7 = 2.9 with a 90% confidence level interval (2.7; 3.2). Figure 5 shows the KS 
significance level Sks vs 7- The curve is very well behaved. This fit requires N s h e ii = 7 ho- 
mogeneous density shells (see eq. 15). We show the density function p(V) for the 7 shells in 
Figure 6; the density rises steeply at large V (last two shells). If we remove these shells from 
the analysis, the results do not change significantly. With 7 = 2.9 our model is consistent 
with both the observed cumulative distribution of N mem and the radial velocity distribution 
with significance levels of SKs(N mem ) — 0-73 and Sk s (V) = 0.62 respectively. 

6. The distribution of velocity dispersions and the group mass function 

To compare with models we derive the velocity dispersion distribution and the mass 
function. Diaferio et al. (1999) make an exhaustive comparison of groups catalogs derived 
from the rCDM and ACDM simulations by Kauffmann et al. (1999) with catalogs derived 
from CfA2N, subset of the UZC. They use the same FOFA we employ and vary the linking 
parameters over a wide range. They also test the sensitivity of their results to variations 
in the luminosity function and to the presence of structures like the Great Wall. They 
demonstrate that for the linking parameters we use here, the catalog derived from the data 
agree reasonably well with the catalogs derived from A CDM provided that the luminosity 
density is the same in both cases. They show that the normalization of the distribution 
of velocity dispersions (and of the mass function) is much more sensitive to the luminosity 
density than to variation in the linking parameters in the group finding algorithm over a 
reasonable range. They emphasize that the Great Wall plays an important role in biasing 
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the distribution of velocity dispersions and the mass function toward higher values relative 
to the simulations. We see a similar effect and comment further below. 

In the previous section we obtain the best fit value 7 = 2.9 (90% c.l. (2.7;3.2)) for the 
exponent of the power-law probability density function of groups, p(N). We conclude that, 
at high confidence level, the density of galaxy systems with true velocity dispersion larger 
than <jt is: 



This distribution applies over the range 100 kms -1 < a < 750 kms -1 . The volume 
covered by the group catalog is too small to include many rich clusters. Within the UZCGG 
there are 20 (robust) groups visible throughout the entire UZCGG volume. These groups 
have N > N r = 840 corresponding to a true velocity dispersion a r = 750 kms" 1 (see eq. 1), 
the actual upper limit of the range of N (and a? ) over which we determine p(N) . For 
the high-<r T systems our analysis does not provide an estimate of the slope of p(N) , but it 
does provides an estimate of the abundance of these systems: [i(N r ) = 1.27 x 10~ 5 /i 3 Mpc -3 . 
We compute the Poisson uncertainty in /z &s(^ r ) in the number of observed groups Nq: 
S/J>obs(N r ) = 0.25 x 10~ 5 /i 3 Mpc~ 3 and conclude that our prediction with the observed value 
Hobs(N r ) to within 0.68 x 5/i b s (N r ). 

We also compare our determination of the abundance of massive systems with previous 
surveys. Direct comparison is possible with Mazure et al. (1996); Fadda et al. (1996); 
Zabludoff et al. (1993). These authors give n(> ax) for clusters. Their system abundances 
at a = 750 kms -1 , are 0.2, 0.6, and 0.6 x 10~ 5 /i 3 Mpc~ 3 respectively. Our system abundance 
is n > 750) = 1.3 x 10~ 5 /i 3 Mpc" 3 . Internal errors are of the order of few tenths in units 
of 10~ 5 . These errors reflect the size of the samples and may be less important than the 
systematic errors introduced by various selection effects. To examine the discrepancy between 
our estimate of n{> 750) and previous ones, we summarize the characteristics of the different 
samples. 

The ENACS clusters (Mazure et al. 1996) lead to the lowest abundance of systems. 
ENACS is a sample of Abell clusters selected according to a richness criterion. Abell clusters 
are an incomplete set of systems (see e. g. Gal et al. 2003, DPOSS2). More importantly, se- 
lection according to Abell richness produces a biased sample of velocity dispersions. Mazure 
et al. (1996) are aware of this problem and discuss it in detail. The velocity dispersions 
of clusters selected according to richness are biased high because a increases with richness 
and because there is a broad scatter around the mean relation. Mazure et al. (1996) use 
reasonable arguments to identify a threshold, a = 800 kms -1 , above which their cluster cat- 



n(> a T ) = (1.27 ± 0.21) x 10~ 5 /i 3 Mpc 



-3 




(23) 
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alog is unbiased. In the end, however, unbiasedness remains an assumption. The abundance 
derived by Fadda et al. (1996) is a factor of three larger than Mazure et al. (1996). Fadda 
et al. (1996) add poorer systems to the ENACS sample. In order to have a a distribution 
that takes into account poor systems, Fadda et al. (1996) scale their richness distribution 
to mimic that of the Edinburgh-Durham Cluster Catalog (EDCC: Lumsden et. al. 1992). 
EDCC is probably a more complete catalog than Abell-ACO, especially at the low richness 
end. Fadda et al. (1996) perform 10,000 random samplings of the velocity dispersions 
of their systems, constraining the richness distribution of the bootstrapped sample to be 
the same as that of EDCC systems. This analysis leads to the same cluster abundance as 
Zabludoff et al. (1993), who use a combined sample of Abell clusters and the thirty densest 
groups selected within the first two slices of the CfA2N redshift survey (Ramella, Geller & 
Huchra 1989). 

We select systems in redshift space from a complete magnitude limited galaxy catalog 
similar to but more extensive than the one used by Zabludoff et al. (1993). In contrast 
with Zabludoff et al. (1993), we include relatively poor systems that would enter neither a 
volume limited sample nor a 2D-selected sample. In fact, because of the large scatter around 
the mean relation between a and richness, several of these relatively poor systems have 
velocity dispersions a > 750). It is therefore not surprising that we find a higher abundance 
of systems above this threshold. 

Other optical and x-ray studies have returned lower system abundances. Bahcall & Cen 
(1993), for example, argue that the abundance of clusters with velocity dispersion larger 
than 750 kms" 1 is 0.2 x 10~ 5 /i 3 Mpc~ 3 , much lower than our estimate. They use a richness 
limited sample of Abell clusters and obtain a scaling relation M oc a 2 significantly different 
from ours (see Table 1). 

The abundance of clusters can also be estimated from x-ray selected samples. The com- 
parison of our abundance of systems with that of x-ray clusters requires a relation between 
x-ray temperature and velocity dispersion and/or mass. The scatter around this relation is 
the source of considerable uncertainty. Reiprich & Bohringer (2002) have carefully analyzed 
the abundance of clusters derived from x-ray data. In the range 2.55 x 10 14 M < M < 
4.4xl0 14 M the cluster abundance is 1.4 xlO" 6 h 3 Mpc~ 3 < n(> M) < 7.0xl0~ 6 h 3 Mpc~ 3 . 
They conclude that the abundance of x-ray clusters is generally lower by a factor of 1.2 — 6.2 
than that of optically selected clusters (Girardi et al. 1998). The origin of the discrepancy 
may be either intrinsic -optical and X-ray clusters may belong to different populations- or 
may result from observational biases affecting the mass estimates. 

We may overestimate n(> 750) somewhat because of (1) projection effects in redshift 
space, (2) large-scale inhomogeneities in the galaxy distribution, and (3) systematic errors 
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in catalog magnitude limits. N-body and geometrical simulations indicate that 20% of the 
5-member groups we study are accidental superpositions. X-ray observations (Mahdavi et 
al. 2000) indicate that at least 40% of the groups we study are true physical systems. 
It is difficult to correct our n(> (Jt) for the presence of false groups since we know their 
abundance nor their velocity dispersion distribution. Whatever the correction, our n(> 750) 
would decrease to come into closer agreement with the estimates of Zabludoff et al. (1993) 
and Fadda et al. (1996). 

Another possible effect driving our abundance toward a high value is the possible over- 
luminosity of the UZC catalog from which we derive our group catalog (Ramella, Geller & 
Pisani 2002). The Zwicky magnitude system may be deeper than implied by its formal 
magnitude limit. In this case, we detect more systems than we should within the surveyed 
volume of the Universe. The density of systems could be too high by as much as a factor two. 
Correction for this effect would bring our estimate of massive systems into close agreement 
with Fadda et al. (1996) and with Zabludoff et al. (1993). 

Finally, it is possible that the region of the UZC is overdense. The abundance of groups 
is proportional to the galaxy density across many different surveys (see e. g. Ramella et al. 
1999). Variations in the density on the scale of the UZC may be as large as a factor of 2, 
but not more. Further, Diaferio et al. (1999) emphasize that the presence of the Great 
Wall in the northern portion of the UZC biases the groups sample toward higher luminosity 
and velocity dispersion relative to a typical regions extracted from a ACDM simulation. The 
groups in the Great Wall are at a larger velocity (~ 8000 kms" 1 ) than the median groups 
redshift (~ 6000 kms -1 ) in a simulated region. These effects may contribute substantially 
to the high normalization of the mass function we derive. 

We next review estimates of the mass function of groups. Girardi & Giuricin (2000) 
measure the mass function of a sample of groups. A direct comparison between their masses 
and ours is not straightforward because their masses are obtained in a model dependent way. 
Taken at face value (see Fig. 9), their mass function is not significantly different from ours. 
Our estimate, however, extends toward lower masses by almost an order of magnitude. 

Heinamaki et al. (2003) estimate the mass function of the Las Campanas Redshift 
Survey groups (Tucker et al. 2000). The amplitude of their mass function is very low (less 
than 0.1 x 10~ 5 /i 3 Mpc~ 3 for masses roughly corresponding to a = 750 kms" 1 . Heinamaki 
et al. (2003) compare their mass function with the results of Girardi & Giuricin (2000) and 
with N-body simulations (see Fig. 7 in Heinamaki et al. 2003). The comparison demonstrates 
the generally low abundance of systems and the marked flattening of the LCRS group mass 
function with respect to the simulated mass function at the low mass end. We obtain a 
similar flattening when we do not take the large-scale variation of group number density into 
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account. Low mass groups are undersampled and a V/V max weighting scheme does not fully 
recover them. 

Martinez et al. (2002) analyze a 2dF sample of groups (Merchan & Zandivarez 2002). 
In this the authors recognize, the flattening of the mass function at low mass results 

from the low number cut-off imposed on their group identification procedure. Martinez et 
al. (2002) compare their mass function with N-body simulations (Jenkins et al. 2001) and 
claim good agreement. Because their best fit agrees with the ACDM model, their result is 
marginally consistent (~ 95%) with ours. 

In Figure 7 we plot n(> ctt) and its 95% confidence level corridor. There we also plot 
the velocity dispersion distribution obtained by Zabludoff et al. (1993) (long dash). As 
expected, the Zabludoff et al. (1993) volume limited sample does not include low velocity 
dispersions systems. 

By using the relations between richness and mass (see Eq. 6 and Table 1), and between 
richness and velocity dispersion (eq. 1) we can convert n(> a?) into an estimate of the mass 
function. We compare this determination with the simulations of Jenkins et al. (2001, 
Virgo Consortium) who derive predictions for the mass function of systems of galaxies over 
our observed mass range. Because the theoretical mass functions predicted by Jenkins et 
al. (2001) apply to masses of critical overdensities, (5p c /p) rn , we scale their masses to 
match our group number overdensity threshold (5p c / ' p)^ = 80. We perform the scaling as in 
Borgani et al. (1999) for both of the cosmological models in Jenkins et al. (2001), ACDM 
{Q m = 0.3, Q A = 0.7, a 8 = 0.9) and rCDM with {Q m = l,Q A = 0, a s = 0.6). 

Figure 8 compares our power law determination of the mass function with the theoretical 
results of Jenkins et al. (2001). We plot the cumulative mass function n(> M) versus the 
mass M. The solid line represents our power law estimate, the dotted lines are 95% confidence 
levels. The rCDM mass functions lies within the 95% confidence level over the mass range. 
The mass function predicted by the ACDM cosmology is only marginally consistent with our 
observations. If we correct the mass function normalization for the presence of unphysical 
groups, the agreement with the ACDM model improves. 

The agreement of our observations with the theoretical mass function is remarkable. 
The use of our weighting procedure is critical in obtaining this result. In fact, if we use the 
standard approach and weight each group by its maximum accessible volume, V/V ma:r , we 
obtain the thin long-dashed curve in Figure 9. This curve shows the characteristic bending 
to a shallower slope at low masses. As discussed above, this shallower slope appears in 
other observed mass functions (e.g.: Martinez et al. 2002; Girardi & Giuricin 2000). Our 
procedure extends the determination of the mass function to masses almost one order of 
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magnitude lower than previous estimates. 

7. Consistency tests and scaling relations 

The mass function we derive depends on the estimate of the two functions p(V) and 
p(N) . We now show that, starting from these two functions, we can reproduce the observed 
distributions of all the main physical parameters of groups, i.e. V, N mem , L mem , a v , and 
virial mass, M vir . 

We proceed in the following way: we randomly sample p(V) and generate a group with 
N members brighter than Lu m . We estimate the number of visible members, N mem , as 
Nmem = N ■ v(V), according to 2. We discard the group if N mem is less than 5. If we 
have at least 5 members, we estimate L me m by randomly sampling the Schechter (1976) LF 
Nmem times. Finally, we estimate a v by randomly sampling a gaussian with dispersion Ot , 
where we compute ot from the best fit model equation 1. We could estimate the virial mass 
by using the true velocity dispersion ot '■ Mt = ?>o\R V i r j G. However, in the usual analysis 
of an observed group catalog, the value of ot is not available and the standard estimate of 
the virial mass is M vir = 3o-^R vir /G. We simulate the observationally derived value of M vir 
in the following way. We compute the scaling relations between, for example, N and M T by 
using the bisector of the two straight lines we obtain first by least squares of log(M T ) versus 
log(iV) and second by least squares of log(iV) versus log (Mr). We compute the simulated 
value of M vir as M vir = M T {(y v / a T ) 2 . 

We stop the previous Monte Carlo simulations when we reach a total of 1000 observable 
groups. We compare the simulated to the observed distributions and compute the signifi- 
cance of the agreement between the two distributions with a KS-test. We repeat the whole 
procedure 10 times and take the median KS significance level as a measure of the agreement 
between the distributions. 

We obtain the following median values of S KS - 0.43, 0.76, 0.29, 0.15, and 0.34 for 
V, Nmem , L me m, &v , and M vir respectively. Figure 10 shows the distributions of the 
physical parameters of 1000 simulated groups (thin line) together with the corresponding 
distributions for the observed groups (thick line). The agreement between the observations 
and the predictions of the simulations indicates that our model is a satisfactory representation 
of the data. 

To test the sensitivity of the exponent 7 in the p(N) model, to the parameters of the N 
vs. <jt law, we use the most extreme values of the parameters a and iVo in Eq. 1 within the 
99% confidence level contour (see Figure 2a). With these parameter values it is impossible to 
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recover both the observed distributions of the group members N mem and velocity dispersion 
a v . 

Table 1 lists the scaling relations we obtain. We determine the scaling relations as the 
angular coefficient and intercept of the bisector of the two straight line fits X versus Y and 
Y versus X, with X and Y any two related quantities in Table 1. We use the results of the 
two fits to characterize the uncertainty in the scaling relations. 

We actually use the relation between mass, M, and richness, N, (Eq. 6) or, equivalently, 
between M and the velocity dispersion, a, to translate n(> cr T ) into a mass function, n(> 
M). We also use this relation to recover the observed distribution of the virial mass M vir . 

Once we establish the relation between M and N, we also have a relation between M 
and luminosity, L. As expected, the slope of log(M) vs. log(iV) and the slope of log(M) vs. 
log(L) are indistinguishable. Girardi et al. (2002) analyze the relation between M and L 
for a different sample of groups and clusters. They obtain M oc £ L34±0 03 ; within about one 
sigma of our relation. 

N-body simulations and semianalytic modeling of galaxy formation performed by Kauff- 
mann et al. (1999) and Benson et al. (2000) gave a reliable prediction of the M/L ratio for 
galaxy systems and also provide an estimate of the dependence of mass-to-light ratio, M/L, 
on L (or M). It is generally found that the ML ratio increases from poor to rich systems, 
and eventually flattens on large scales. 

Yee & Ellingson (2002) measure the relation between mass and richness (B cg ) within 
the CNOC1 sample of clusters of galaxies. These authors find M oc fi^ 64±0 ' 28 , in very good 
agreement with our result. Because the intercept of this relation depends on the magnitude 
completeness limit, we can compare only slopes and not intercepts. 

As a consistency check, we also compute the relation between log(L^) and log(<T T ). 
Reassuringly, the slope of this scaling law is the same as our best fit of the power law model 
in Eq. 1, iV vs. cr T . 

Finally, for the sake of completeness, in Table 1 we list the relations between log(M) 
vs. log(f7r), log(iV) vs. log(R vir ) and \og(a T ) vs. \og(R vir ). Yee & Ellingson (2002) fit the 
scaling between the virial radius and richness and find B cg oc i? 21±a7 , in very good agreement 
with our result. Our scaling relations significantly extend the data range spanned by cluster 
analyses. These scaling relations are remarkably consistent with those derived for clusters. 

We also explore different forms for p(N). For example, we consider a p(N) that gives an 
exponential law for /(or), i.e. the law used by Zabludoff et al. (1993). This form for p(N) 
fails to reproduce the space density of robust groups /i f, s (> N r ) and leads to a discordant 
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distribution of L mem . Another model we test is the PS model. The problem with this model 
(Equation 7) is that it has two free parameters, n e ff and N+: our procedure can constrain 
only one parameter (see section 6). However, we can test the relation between n e ff and if 
p(N) can be approximated by a power-law (see Equation 10). We compare the PS model for 
p(N) with our data using a KS test. In Figure 4 we plot the 10%, 50%, and 99% confidence 
level contours in the n e // — plane. In the same plane we also plot as a function of 
n e ff according to Equation 10 where we assume < N >= 300, close to the average value for 
our sample. For 7 we use our best fit value, 7 = 2.9. Figure 4 clearly shows good agreement 
between the power-law approximation and the PS model (note that the rejection region is 
external to the outermost contour (c.l. 99%)). 

Sato et al. (2000) determine n eff = -1.2 ± 0.3 for masses in the range 10 12 - -1O 15 M 
based on the analysis of an ASCA x-ray cluster sample. Zandivarez, Abadi & Lambas 
(2001) find similar results for these scales. If we use the value n e // = —1.2 in Equation 7, 
we obtain a PS model with N± ~ 280 corresponding to M* = 1.6 x 10 14 h^ 1 M Q ( 95% 
confidence interval goes from 6.5 x 10 13 h^ 1 M Q to 3.6 x 10 14 h^ 1 M Q ). This value of 
M* is in reasonable agreement with the determinations of Bahcall & Cen (1993), namely 
M* = (1.8 ± 0.3) x 10 14 h' 1 M and Girardi et al. (1998): 2.61°;® x 10 14 h' 1 M . This 
result indicates that our analysis of low-mass systems of galaxies leads to an estimate of the 
mass function parameters consistent with the results determined in the higher mass range 
of galaxy clusters. 



8. Summary 

We measure the probability density function for the velocity dispersion of groups in 
the UZC galaxy catalog. Our method a) does not require that the group spatial density is 
constant, b) predicts the observed distributions of all the fundamental quantities of groups, 
i.e radial velocity, number of members, luminosity, velocity dispersion, and virial mass, and 
c) includes low luminosity systems that would be neglected with the usual analysis of a 
volume limited sample. The best fit for the slope of n(> ax) is —3.4 (—2.6; —4.7) over the 
range 100 kms -1 < a < 750 kms -1 . Over this range, the Press-Schechter model is well- 
approximated by a simple power law model. Our model predictions agree quite well with 
the predictions of Jenkins et al. (2001) for the CDM models. 

Our weighting procedure is critical in the mass function determination. If we use the 
standard approach and weight each group by its maximum accessible volume, V/Vmax, we 
obtain the characteristic bending of the mass function to a shallower slope at low masses. 
This shallower slope appears in other observed mass functions (e. g. Martinez et al. 2002; 
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Girardi & Giuricin 2000). Our procedure extends the determination of the mass function 
to masses almost one order of magnitude lower than previous estimates. 

For massive systems, the density of groups with true velocity dispersion larger than 
750 kms -1 is (1.27 ± 0.21) x 10 5 /i 3 Mpc" 3 . Our determination is larger by a factor 
of two than previous determinations. We examine this discrepancy and identify several 
possible sources for this discrepancy: a) incompleteness of cluster catalogs used in previous 
works, b) presence of false groups in our catalog produced by projection effects in redshift- 
space, c) systematic errors of the Zwicky magnitude system (may be deeper than its nominal 
magnitude limit would imply), d) real overdensity of the region of the UZC Corrections for 
any of these effects would bring the abundance determinations into closer agreement. 

The mass function we derive depends on the estimate of the two functions p(V) and 
p(N) . Starting from these two functions, we reproduce the observed distributions of all 
the main physical parameters of groups, i.e. V, N mem , L mem , a v , and virial mass, M vir . 
We also obtain scaling relations between these physical parameters. Our scaling relations 
significantly extend the data range spanned by cluster analyses and are remarkably consistent 
with those derived for clusters. 
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Fig. 1. — The relation between the true velocity dispersion o T and observed velocity disper- 
sion a v for a Monte Carlo simulation. The solid line represents the identity. 
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Fig. 2a. — The contour levels of the model N versus ot (eq. 1) for the UZCGG groups 
with 5 or more members. The outer thick contour traces the 99% confidence region, the 
middle contour traces the 90% confidence level, and the inner thin contour traces the 50% 
confidence level. 
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Fig. 2b. — Same as Figure 2a for x-ray emitting groups 
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Fig. 3. — The relation between the group richness N and the velocity dispersion a v for 
the x-ray emitting groups. The thin line results from the bisector of the two least squares 
straight lines (one by fitting iV vs. a v and the other by fitting a v vs. N), the thick line 
shows that our Monte Carlo approach has less sensitivity to low a v groups. 
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Fig. 4. — The PS model: the contours show the KS test significance level SKs{N mem ) for 
the comparison between the observed distribution of the groups and the PS model described 
by N* and the effective spectral index n e ff. The outer contours indicate the 99% confidence 
level, the middle contours indicate the 50% confidence level and the inner contours show the 
10% confidence level. The dot-dashed line shows the relation for data distributed according 
to a power law (eq. 10) and with < N >= 300, 7 = 2.9 and k = 1.43. The power-law line 
lies well within the 50% confidence level. The relation between the group richness iV and 
mass is M oc N 1A3 (see Table 1). 
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Fig. 5. — The power law model: the KS test significance level SKs(N mem ) for the comparison 
between the model and observed distribution of the group members. 
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Fig. 6. — [Fig. 6] The density of groups as a function of the radial distance p(V). The 
function is the best fit power law model for p(N) with = 7 intervals. The error bars in 

V are the width of the bins; for p they the Poisson uncertainty 5Ni in the number of groups 
per bin. 
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Fig. 7. — The power law model for the space density of groups with true velocity dispersion 
> a T . The thin curve shows the power law model (Eq. 23), the dotted lines show the 95% 
confidence band for the entire range of ot- The right side vertical bar shows the width of 
the 95% confidence level for the observed density of robust groups. The dashed curves show 
the data of Zabludoff et al. (1993) (± 1 — a poissonian error band). 
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Fig. 8. — The cumulative mass function of groups (thick solid line). Dotted lines are 95% 
confidence levels. The thick short- and long-dashed lines are theoretical predictions (Jenkins 
et al. 2001) for ACDM (Q m = 0.3, fi A = 0.7, a 8 = 0.9) and rCDM with (Q m = 1,Q A = 
0, <7 8 = 0.6) models respectively. 
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Fig. 9. — The cumulative mass function of groups (thick solid line). Dotted lines are 95% 
confidence levels. The dashed line shows the results obtained by Girardi & Giuricin (2000). 
The thin solid curve is the observed mass function we obtain assuming constant large scale 
density of groups. 
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Fig. 10a. — [Figure 10a] The Power Law model for the group pdf: the observed and model 
cumulative distribution of radial velocity V of UZCGG groups (thick curve) and a 1000 
group Monte Carlo simulation of our best fit power law model (thin curve). 




Fig. 10b. — [Figure 10b] The power law model for the group pdf: observed and model 
cumulative distribution of N mem of UZCGG groups (thick curve) and a 1000 group Monte 
Carlo simulation of our best fit power law model (thin curve).] 
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Fig. 10c. — [Figure 10c] The power law model for the group pdf: observed and model 
cumulative distribution of log(L mem ) of UZCGG groups (thick curve) and a 1000 group 
Monte Carlo simulation of our best fit power law model (thin curve).] 
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Fig. lOd. — [Figure lOd] The power law model for the group pdf: observed and model 
cumulative distribution of a v of UZCGG groups (thick curve) and a 1000 group Monte Carlo 
simulation of our best fit power law model (thin curve).] 
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Fig. lOe. — [Figure lOe] The power law model for the group pdf: observed and model 
cumulative distribution of virial mass M„j r of UZCGG groups (thick curve) and a 1000 
group Monte Carlo simulation of our best fit power law model (thin curve).] 
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Table 1. The scaling relations among the main group parameters. 



Relation Slope Intercept 

log(M T ) vs. log(iV) 1.43±g;8t 10.70^;°9 

log(M T ) vs. log(L T ) 1.44±g;8i -2.2ltg;?? 

log(M r ) vs. \og(a T ) 2.58l°;°? 7.48^ 

log(Lr) vs. log(<7 T ) 1.80+°;^ 6.71±°;^ 

log(iv) vs. io g (i^ r ) 2.ool:lii 2.u-§™ 

log(a T ) vs. log(itW) 1.18+^5 2.42^ 



